A theory of jet shapes and cross sections: from hadrons to nuclei 
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For jets, with great power comes great opportunity. The unprecedented center of mass energies 
available at the LHC open new windows on the QGP: we demonstrate that jet shape and jet cross 
section measurements become feasible as a new, differential and accurate test of the underlying 
QCD theory. We present a first step in understanding these shapes and cross sections in heavy 
ion reactions. Our approach allows for detailed simulations of the experimental acceptance/cuts 
that help isolate jets in such high- multiplicity environment. It is demonstrated for the first time 
that the pattern of stimulated gluon emission can be correlated with a variable quenching of the jet 
rates and provide an approximately model-independent approach to determining the characteristics 
of the medium-induced bremsstrahlung spectrum. Surprisingly, in realistic simulations of parton 
propagation through the QGP we find a minimal increase in the mean jet radius even for large jet 
attenuation. Jet broadening is manifest in the tails of the energy distribution away from the jet axis 
and its quantification requires high statistics measurements that will be possible at the LHC. 

PACS numbers: 



I. INTRODUCTION 

Fast partons propagating in a hot/dense nuclear 
medium are expected to lose a large fraction of their 
energy In fact, the stopping power of strongly- 

interacting matter for color-charged particles has, by far, 
the largest experimentally established effect: the atten- 
uation of the cross section for final-state observables of 
large mass/momentum/energy. This jet quenching mech- 
anism has been used to successfully explain the strong 
suppression of the hadron spectra at large transverse 
momentum observed in nucleus-nucleus collisions at the 
Relativistic Heavy Ion Collider (RHIC) 0,0]. There is 
mounting evidence that the quark-gluon plasma (QGP)- 
induced quenching can be disentangled from other nu- 
clear effects even at the much lower Super Proton Syn- 
chrotron (SPS) center of mass energy [J]. To calculate 
parton energy loss in the QGP several theoretical ap- 
proaches have been developed d, d, 0, Hi ■ While tackling 
the same basic problem, they use different assumptions 
for the boundary conditions (initial/final quark or gluon 
virtuality), make different approximations for the parent 
parton and radiative gluon kinematics, and treat differ- 
ently the interaction between the jet-|-gluon system with 
the medium (differentially vs on average). For discussion 
see 0. 

At present, most measurements of hard processes are 
limited to single particles and particle correlations, which 
are only the leading fragments of a jet. There is general 
agreement on the physics that controls inclusive particle 
suppression in the QGP and the experimental methodol- 
ogy of determining Raa {pt ) (or Iaa {pti ,Pt^)) [2] . Thus , 
for any particular combination of radiative/coUisional en- 
ergy loss evaluation and its phenomenological application 
to leading particle quenching the QGP density can be de- 
termined with 20 — 25% "statistical" accuracy [13. An 



inherent limitation of this approach is that while fits to 
the data can, not surprisingly, always be performed [Tlj 
they do not resolve the staggering order of magnitude 
" systematic" discrepancy in the extracted medium prop- 
erties. Furthermore, focusing on quantities that can be 
constrained with little ambiguity from the measured ra- 
pidity entropy density in heavy ion collisions distracts 
from issues such as the approximations that go into the- 
oretical energy loss derivations and their application to 
systems where more often than not these initial assump- 
tions are violated [13] ■ In searching for experimental 
measurements which can pinpoint the framework for en- 
ergy loss calculations that is applicable to heavy ion re- 
actions, complex multi-particle correlations may not be 
optimal. They are very sensitive to non-perturbative ef- 
fects/fragmentation and the modeling effort cannot be 
systematically improved due to the violation of factoriza- 
tion for highly exclusive observables [13. It is, therefore, 
critical to find alternatives that accurately reflect the en- 
ergy flow in strongly-interacting systems and have a more 
direct connection to the underlying quantum chromo- 
dynamics (QCD) theory. 



The intra-jet energy distribution and the related cross 
section for jets in the case of heavy ion reaction closely 
match the criteria outlined above. The high rate of hard 
probes at the LHC and the large-acceptance calorimetry, 
see e.g. [l3|, will enable these accurate measurements. It 
should be noted that proof-of-principle measurements of 
jet cross sections have become possible at RHIC [3], but 
significantly better statistics will be required to quantify 
the QGP effects on jets. In this paper we study the 
magnitude of these modifications in Pb-|-Pb collisions at 
-^i = 5.5 TeV at LHC. We demonstrate that a natural 
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generalization of leading particle suppression to jets, 



pjct / TP^. Tjmax mm 



) = 



d<jPP(ET\R" 



(1) 



dyd^Er 



is sensitive to the nature of the medium-induced energy 
loss. The steepness of the final-state differential spectra 
amplifies the observable effect and the jet radius i?™^'^ 
and the minimum pRrticle/tower energy pj^ min ~ cj^^^^ 
provide, through the evolution of i?-j^^(£'T; i?™'''', w""") 
at any fixed centrality, experimental access to the QGP 
response to quark and gluon propagation. 

In the following discussion of jet shapes and jet 
cross sections in p-fp and A-f A collisions we stay as 
close as possible to an analytic theoretical approach. 
Thus, we are able to unambiguously connect the non- 
perturbative QCD effects, medium properties and the 
induced bremsstrahlung spectrum to experimental jet 
observables. Determination of baseline jet shapes and 
their generalization to finite momentum acceptance cuts 
builds upon the work of Seymour [l^. We refer the 
reader to [Tsl . [l^ for discussion of the complications in 
defining a jet and the related topic of jet-finding algo- 
rithms. To make the discussion simpler, we will assume 
that the complications of the different definitions can be 
subsumed into an Rsep parameter, as described in Ap- 
pendix [BJ Once a jet axis and all of the jet particles / 
calorimeter towers "i" have been identified, the "integral 
jet shape" is defined as: 



*i„t (r-;i?) 



E,(i^T).9(r-(j?jet).) 



(2) 



where r, R are Lorentz-invariant opening angles, Rij — 

\J iVi ~ Vj)^ + {4'i ~ 4>jY 1 ^iid i represents a sum over 
all the particles in this jet. ^int(?';^) is the fraction 
of the total energy of a jet of radius R within a sub- 
cone of radius r. It is automatically normalized so that 
^'i„t(i?; R) = 1. To move from the integrated to the dif- 
ferential jet shape, we define: 



ij{r-R) = 



rf^int(7-;fl) 

dr 



(3) 



This is the angular density of jet energy (remember- 
ing that the appropriate 3D representation would be 
■0"*^(r;i?) = ■^tp{r] R)). Understanding the many-body 
QCD theory behind jet shape calculations will naturally 
lead to understanding the attenuation of jets in reactions 
with heavy nuclei. 

This article is organized as follows: in Section |TT1 we 
outline a calculation of the jet shape in nucleon-nucleon 
(N-l-N) collisions using the framework of perturbative 
QCD. We compare this calculation to existing Tevatron 
data and investigate the jet shapes at LHC energies. A 
brief discussion of final-state QGP-induced radiative en- 
ergy loss in the GLV formalism is given in Section IIIII 
We prove that the cancellation of small-angle near-jet 



axis bremsstrahlung persists to all orders in the correla- 
tion between the multiple scattering centers and provide 
details of its numerical evaluation. The fully differential 
distribution of the energy lost by a hard parton is also 
shown. In Section HVl we present results for the medium- 
modified jet shapes and cross sections as a function of the 
jet cone radius, R, and the experimental pt cut, Wmin, 
and discuss a simple energy sum rule. We demonstrate 
the connection between the characteristic properties of 
the QGP-induced gluon radiation and the variable sup- 
pression, at the same impact parameter h, of jet rates, 
the modulation of the mean jet radius and the enhance- 
ment in the "tails" of the intra-jet energy flow distri- 
bution. A summary and conclusions are presented in 
Section |Vl Appendix [X] shows a calculation of the base- 
line jet cross sections at the LHC and an estimate of the 
accuracy with which these cross sections and jet shapes 
can be measured with nominal first-year integrated lumi- 
nosities in p+p and A+A reactions. In Appendix iBl we 
study the influence of the different perturbative and non- 
perturbative contributions to the jet shape in hadronic 
collisions. Finally, Appendix [Cl contains a discussion of a 
double differential measure of energy flow in jets and its 
connection to particle angular correlation measurements, 
currently conducted at RHIC. 



II. JET SHAPES IN 'ELEMENTARY' P-P 
COLLISIONS 

In the process of advancing perturbative QCD theory 
to many-nucleon systems, the final-state experimental 
observables should be first understood in the simpler p-fp 
reactions. 



A. Theoretical considerations 

1. Leading order results 

In the introduction, we defined the central quantity 
of our study, the differential jet shape for parton "a", 
"fpaii", R)- As in [l^, the starting point of the calculation 
is the leading order parton splitting: a suitable separa- 
tion of physical time scales enables the separation of the 
calculation into production and jet showering. The QCD 
splitting functions Pa^bc{z) give the distribution of the 
large fractional lightcone momenta (or approximately the 
energy fractions) of the fragments relative to the parent 
parton, z and 1 — z respectively. To lowest order, recall- 
ing that Tpair] R) describes the energy fiow oc 2, we can 
write: 



rf^'i„t,a(r;i?) a. 2 '•1-^ 



dr 



(4) 

In Eq. (HI r = (1 — z)p is related to the opening angle p 
between the final-state partons. 
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In "elementary" p+p collisions the inclusion of soft 
particles {zmin ~ 0) in theoretical calculations is not a 
bad approximation. Even in this case, however, there are 
intrinsic limitations on the minimum particle/calorimeter 
tower pt or Et, related, for example, to detector accep- 
tance. In heavy ion reactions, especially for the most 
interesting case of central collisions, there is an enor- 
mous background of soft particles related to the bulk 
QGP properties. Jet studies will likely require minimum 
particle energy > 1 — 2 GeV at RHIC and even more 
stringent cuts at the LHC. Furthermore, control over 
Zmin can provide detailed information about the proper- 
ties of QGP-induced bremsstrahlung. Further kinematic 
constraints on the values of z arise since both the result- 
ing partons must be within an angular distance R of the 
original jet axis, r < R, rz/{l — z) < R. In this case 
they are identified with the jet. If not, they are iden- 
tified as two separate jets. For a cone-based algorithm, 
the relative separation i?sep^ (as opposed to just the dis- 
tance from the original jet axis) is an additional criterion: 
p < RsepR- We find: 



In contrast to the case of e~^+e~ annihilation, hadronic 
scattering is accompanied by copious initial-state radia- 
tion (ISR) that can fall within the jet cone. While the 
contribution of the ISR is small for small values ofr/R, it 
gives an essential contribution at larger angles. A simple 
estimate based on a dipole radiation and the kinemat- 
ics of the hard parton - soft gluon coincidence within a 
cone [l5l|, similar for both quark and gluon jets, yields: 



(9) 



Again, the 'plus-prescription' to account for the r = 
point is not explicitly shown. In Eq. ([5]) C ~ Cp ~ CaI'^- 
The leading order calculation is most appropriate for 
the rare, hard splittings of a very high momentum jet. 
The use of a running coupling improves the numerical 
results, providing larger weight for softer events. We 
employ a running as evaluated at the largest fcy in 
the problem, [i = r(l — Z)Et for the jet splitting and 
fi — (1 — Z)Et for the initial state radiation |15i |. 



Z = max < z, 



Z = max < z 



r + R 



RsepR 



if r < {Rsep - l)R , (5) 



if r > [Rsep - l)R . (6) 



Carrying out the integration in Eq. (|4|) we arrive at the 
LO jet shape functions for quarks and gluons: 



CpOls 2 / 1 - Zmm 

Mr) - ^7 21°g^^ 



-^[(l-Z)2-zL„]) , 



(7) 



V'gW = -T^Z 1 2 log- 



27r r 



Z 



2z^ ■ - -2^ ■ + -z'^ ■ 



-z^ ■ +z^ ■ 

o min ' min 



(8) 



In the Zjnin limit Eqs. ([7]) and ([5]) reduce and we 
recover the previously known results |15j. There is an 
implicit 'plus-prescription' in these results when calcu- 
lating moments of physical quantities, as we have not 
considered the virtual corrections in the forward direc- 
tion. Hence, the result is not applicable for r — and 
does not have the correct normalization when integrated. 
However, the shape is reflective of final-state parton split- 
ting and, when needed, for the leading order calculation 
one may apply a cutoff for small r and normalize via 
first-bin subtraction. 



2. Resummation - all orders and multiple emission 

As r ^ 0, in the coUinear limit of parton splitting, the 
leading order contributions to the jet shape diverge, see 
Eqs. (O, ^ and In fact, all orders in the perturba- 
tive expansion diverge, including powers of logr in the 
form Q;"log^"^^r. With plentiful parton showering, it 
becomes increasingly less likely that any particular quark 
and gluon will be coincident with the jet axis. Quanti- 
tatively, this is described by a Sudakov form factor. The 
energy density at small angles is dominated by the hard 
parton in the splitting. If there is a splitting that leaves 
the hard parton at an angle ri , a subsequent splitting at 
r2 < Ti will not contribute to the energy density at r2. 
Multiple independent splitting follows a Poisson distri- 
bution, hence, the probability of energy flow at an angle 
less than r is exponentially suppressed by the integrated 
probability at angles greater than r, i.e.: 

P(< r) = exp(-Pi(> r)) (10) 
= exp (- / dr>coii(r') I • (H) 



This only applies for soft emissions which do not take 
away (much) momentum, i.e. at leading log accuracy. 
Improvement can be obtained at modified leading log ac- 
curacy (MLLA), when the running of the coupling con- 
stant is included in P\. We don't take other (e.g. recoil or 
kinematic constraints) effects in evaluating the Sudakov 
form factor in the soft coUinear approximation. 

The resummed V'rcsuin('') = 'jp-^i'^ carry 
out the integration in Eq. including the running 

asirEx), to obtain modified leading logarithmic accu- 
racy (MLLA). Note that as(/z) = l/(2/3o log ^^-^), with 



47r/3o = bo 



^-Ca - ^TnNf. First, take the small r 
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limit in Eqs. (O, ([8]) and ([9]), keeping terms oc 1/r and 
oc l/rlog(l/r). Based on Z — max(zmi„, r/i?), in this 
limit we have two kinematic domains. For r > ZminR the 
results are similar to the known case of no acceptance 
cut-off and reduce to the known results if Zmin = 0: 



R 



Pq{r > z,ninR) = exp 2Cf log — /i 2/3oas log — 



R 



-Cp — CR? — C^{Zmin) 



R 



X /2 2/3oas log — 



(12) 



R 



Pg{r > ZrainR) = exp 2C A log — /i 2/3oas log — 



R 



-6o - CR^ - C^{Zratn) 



X /2 2/3oas log 



i? 



(13) 



It is useful to employ the same notation as in ^15i] and fa- 
cilitate the comparison to the case of no kinematic cuts: 
/i(x) = log(l - a:)/(27r/3o), and ^(x) = (1 - log(l - 
x) / x) / {2t: (3q) . The Zmin- dependent corrections are iso- 
lated as follows: 



C-n ZfriinR-i Zmin 

) = 2Cf l0g(l - Zrmn) 

-^C.z' 



(14) 



^min-^-) 



2Ca log(l - z„ 



+Ca 2z, 



2 mm 
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(15) 



When r < ZmmR the integration in Eq. (|lip has to be 
split in two regions: r' G {r,ZminR) and r' G {zminR, R)- 
The second integral is trivially obtained from the case 
that was just considered above, Eqs. (I12p and with 
the substitution r = ZminR- When combined with the 
first integral, it yields: 

Rqir ^ ZjyiinR) Rqir ^ Z^inR^^ — ZminR) 



X exp 



x/2 2/3oas log- 



(16) 



Rgi^ ^ ZminR) Rgi^ ^ ZminRi"^ — ZmijiR) 



X exp 



-&0 - C>{Zmin) 



x/2 2/3oa3 log- 



(17) 



Here, we denote by = asizminREr), as opposed to 



as = as{RET), and: 



C"^ {r < ZmmR; Zmin) = 2C_F loi ' ""^ 



(18) 

C;^(»' < ^min^?; Zmin) = 2Ca log ""^ 



~Tj^Nf (^Zmin ^ "^^min + g^min^ ■ (19) 

Note that the Sudakov form factors evaluated here will 
regulate any coUinear divergence present in V'coU- 



3. Power corrections - an estimate of non-perturbative 
effects 

Inclusion of the running coupling constant under the 
momentum transfer integrals yields contributions from 
regions in which Q ^ ^qcd or lower, i.e. there is a 
fundamental non-perturbative contribution to all of the 
integrals. An estimate of these power correction effects 
with finite acceptance gives the following result: 



i^pc{r) 



2Cr 2 Qo 



ao'(<9o, k„ 



2tt r tEt 



wo 



+- 



2Cr 2 k„ 



27r r tEt 



as{p) 



+2|3oasi^l)^ 1 + log 



kr, 



(20) 



where Cr = Cf,Ca for quarks or gluons, respectively. 
In Eq. ((20| kmin — Zmin^Ex and /i is the renormaliza- 
tion scale. The term oc cib'(Qo, fcmm) depends on the 
parametrized non-perturbative contribution, defined as: 



1 f'^" 

aa'{Qo,km2n) = -pr- dkasik) , (21) 

with Qo representing the non-perturbative scale. In our 
numerical calculation we use 

do' {2 GeV, 0) = 0.52 , ao'(3 GeV, 0) = 0.42 (22) 

from Ref. [3, [l3| and a parametrization of the strong 
coupling constant at small momentum transfer given in 
Ref. 119|. The terms oc as(M)i ^^(/i), come from subtract- 
ing the perturbative component in the non-perturbative 
region [18]. Finally, the term oc kmin results from the 
introduction of finite acceptance, Zmin- 
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A similar expression is derived for initial-state radia- 
tion: 



'4'i,pc{r) 



— 2r— Uo (Qo,fc„ 



-2/3oas(A')Ml + log7^ 



as{^JL) 



Qo 



2tt 



+2fioas{pLf l + lo, 



, (23) 



where = ZminEr, and C ~ Cp ~ Ca/2. At lower 
jet energies the power corrections are, in fact, sizable 
even at large r. This suggests that if there is deviation 
between the theoretical results and the experimental data 
it may be largely due to incomplete consideration of non- 
perturbative effects. We stress that the generalization of 
power corrections to finite acceptance implies that these 
should be taken into account only for ZminfET < Qo or 
ZminET < Qo for fiual-statc and initial-state radiation, 
respectively. 



4- Total contribution to the jet shape 

As indicated before, the resummed jet shape at small 
r/R is evaluated as i/'rcsum(r) = -^P{r) = ^coii{r)P{r). 
Taking all contributions to the jet shape and ensuring 
that there is no double counting at small r/R to 0{as) 
we find: 

i'ir) = V-couW (P(r) - 1) + VLo(r) +V'i,Lo(0 

+^PcW +V'«,PcW , (24) 



On the right-hand-side of Eq. (|24p the first term comes 
from Sudakov resummation with subtraction of the lead- 
ing 1/r, (1/r) log(l/r) contribution at small r/R to avoid 
double counting with the fixed order component of the 
differential jet shape. The second and third terms repre- 
sent the leading-order contributions in the final-state and 
the initial-state. The last two terms represent the effect 
of power corrections. In a full calculation the relative 
quark and gluon fractions fq + fg = 1 are also needed: 
i^{r;ET) = MET,V^)Mr;ET) + fg{ET,V~s)Mr;ET). 
These fractions are calculated in Appendix [X] alongside 
the demonstration of the feasibility of jet cross section 
and differential jet shape measurements at ^/s — 5.5 TeV. 

If the resummed part completely dominates the area 
under ip{r) in Eq. (|24p . i.e. the power corrections and the 
fixed order result affect only the large r/R "tails", the 
theoretically calculated differential jet shape is properly 
normalized. In reality, this is not the case and the first 
correction, 0{al), arises from '0con('')(-P(< r) — 1) = 
'0coii('')(l + C'ii;as(- • • )-!-• — 1). It will be larger for gluon 
jets when compared to quark jets, Ca vs Cp^ and for 
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FIG. 1: (Color online) Comparison of numerical results from 
our theoretical calculation to experimental data on differential 
jet shapes at = 1960 GeV by CDF II Insert shows 

the Et dependence of Rsep- 



lower transverse energies. The normalization can them 
be ensured via 

iP{r) ^{r) + Norm x Vcon(?-) ln(P(< r)) , (25) 

where Norm is determined numerically. We stress that 
to achieve a robust theoretical description of the differ- 
ential jet shape all contributions to ^/'(r) from Eq. 
should be included. In Appendix |B] we elucidate their 
relative strength using numerical examples. We also in- 
vestigate the dependence of the shape on Rsep and the 
non-perturbative scale Qo- 



B. Comparison to the Tevatron data 

In Fig. [T] we show comparison of the theoretical model 
for the jet shape, Eq. p4|) . to the experimental measure- 
ments in p + p collisions at ^/s = 1960 GeV at Fermilab 
from Run II (CDF II) 2"^. Our numerical results include 
all contributions from leading order, resummation and 
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power corrections with Qo = 2 GeV. The insert shows 
the variation of the parameter Rgep with the transverse 
energy of the jet. At high jet Et our theoretical model 
gives very good descriptions of the large r/R experimen- 
tal data with i?sep = 1.3-1.4. For Et = 45-55 GeV the 
largest meaningful value Rsep — 2 can describe the data 
fairly well, except at very small r/R region. Extended 
discussion of the various contributions to the differential 
jet shape is given in Appendix [B) 

We note that for r/i? <C 1 and a large gluon jet frac- 
tion in conjunction with moderate Et < 50 GeV there is 
still deviation between the data and the theory, e.g. the 
top panel of Fig. [1] This is likely related to the need for 
significant corrections, Eq. (j25p . to ensure the proper nor- 
malization of 'ip{r/R). Such corrections, in turn, point to 
NLO effects, a possible breakdown of our soft coUinear 
jet splitting approximation for the Sudakov resumma- 
tion and non-perturbative effects. Note that even Monte 
Carlo event generators have to be tuned to describe this 
data [13] . For the purpose of our manuscript the deficien- 
cies in this specific part of phase space are not essential 
since, as we will see in Section llVl the experimental signa- 
tures of jet propagation in the QGP are most pronounced 
in the complementary r/R 1 domain. One simply 
has to keep in mind that the description oi r/R < 0.25 
^{r / R) in the vacuum allows for further theoretical im- 
provement. 



C. Predictions for the LHC 

We employ the theoretical model that describes the 
CDF II data and apply the same transverse energy- 
dependent Rsep parameter to obtain predictions for the 
LHC at ^/s = 5.5 TeV. The emphasis here is to produce 
a baseline in p -\- p reactions for comparison to the full 
in-medium jet shape in Ph + Pb collisions. The essential 
diff'erence in going from the Tevatron to the LHC is in the 
production of hard jets. At the higher collision energy we 
observe a greater contribution from gluon jets relative to 
quark jets, e.g. Fig. [12] in Appendix [Xj Therefore, for 
the same Et^ jets at the LHC are expected to be slightly 
wider than at the Tevatron. 

Figure [2] shows our numerical results for the jet shape 
for four different energies Et = 50, 100, 250 , 500 GeV 
and two cone radii R — 0.7, 0.4 in p + p collisions at 
y/s = 5.5 TeV at LHC. An interesting observation is 
that, when plotted against the relative opening angle 
r/i?, these shapes are self-similar, i.e. approximately in- 
dependent of the absolute cone radius R. One of the main 
theoretical developments in this paper is the analytic ap- 
proach to studying finite detector acceptance effects or 
experimentally imposed low momentum cuts. In Fig. [2] 
this is illustrated via the selection of Zmm = Pt mm /Et — 
0.2, 0.1, 0.04 and 0.02 {pT min = w™" = 10 GeV). Elim- 
inating the soft partons naturally leads to a narrower 
branching pattern. However, for this effect to be read- 
ily observable 10 — 20% of the jet energy, going into soft 
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FIG. 2: (Color online) Numerical results for the differential jet 
shapes in p-l-p collisions at -^s = 5.5 TeV at the LHC. Solid 
lines represent jet shapes with R = 0.7, a;™'° = GeV, dashed 
lines stand for jet shapes with R = 0.4, w™'" = GeV, and 
dashed-dotted lines are for jet shapes with R = 0.7, oj™'" = 
10 GeV. The inserts show integrated jet shapes 'I'int(r; _R). 



particles, must be missed. Thus, even with pT mm ~few 
GeV cuts in Pb-|-Pb collisions at the LHC aimed at reduc- 
ing or eliminating the background of bulk QGP particles 
that accidentally fall within the jet cone, the alteration of 
ip{r/R) is expected to be small. We also studied the in- 
tegral jet shape ^'int(r/i?), shown in the inserts of Fig. O 
as a tool for identifying kinematic and dynamic effects 
on jets [l3l- Only when large differences exist between 
two il}{r/R) for r/R < 0.4 these will be refiected in the 
integral jet shape. If the differences are pronounced in 
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FIG. 3: (Color online) 3D plot of the diflterential jet shapes at 
three different jet energies Et ~ 20 GeV (top panel), Et = 
100 GeV (middle panel), and Et = 500 GeV (bottom panel) 
with R = 0.7, = GeV in p+p collisions with ^/s = 

5500 GeV at the LIfC. From low jet energy to high jet energy, 
jet shape becomes much steeper. 



the r/R > 0.4 region, as is typically the case for heavy 
ion reactions, ^'int(?'/^) will be practically insensitive to 



the QGP effects on jet propagation. 

It is important to also note in Fig.[2]that there is a dra- 
matic change in the differential jet shape in going from 
small to large transverse energies even in p + p reactions. 
This is best demonstrated on a 3D-plot, where the vol- 
ume of the jet cone is normalized to unity. Examples of 
'^tildr'^ for Et = 20, 100, 500 GeV are given in Fig. H 
We have used R = 0.7 and w™™ = 0. It is obvious 
that development of detailed theoretical models and their 
validation against experimental data in nucleon-nucleon 
collisions are necessary before any credible conclusions 
about the modification of the QCD jets in the QGP 
medium can be drawn. 



III. MEDIUM-INDUCED CONTRIBUTION TO 
THE JET SHAPE 

The principal medium-induced contribution to a jet 
shape comes from the radiation pattern of the fast quark 
or gluon, stimulated by their propagation and interaction 
in the QGP. There is a simple heuristic argument which 
allows one to understand how interference and coherence 
effects in QCD amplify the difference between the en- 
ergy distribution in a vacuum jet and the in-medium jet 
shape [2l[. Any destructive effect on the integral av- 
erage parton energy loss Ai?''"'', such as the Landau- 
Pomeranchuk-Migdal effect, can be traced at a differen- 
tial level to the attenuation or full cancellation of the 
collinear, uj, gluon bremsstrahlung: 



A prad 

'-^^LPM suppressed 



dI3 

duj 



[lU ^ -E)lPM suppressed 



d/9 



dojd^kT 



LPM suppressed ; 



(26) 



and we indicate the parts of phase space where the mod- 
ification of the incoherent dl^ /duid^kr is most effective. 

Indeed, detailed derivation of the coherent inelastic 
parton scattering regimes in QCD was given in In all 
cases, the origin of the LPM suppression can be traced 
to the cancellation of the collinear bremsstrahlung. 
The destructive quantum interference is most promi- 
nent for final-state radiation, where the large-angle gluon 
bremsstrahlung was originally discussed in Ref. 22] to 
first order in opacity. Even though to carry out realistic 
simulations to higher orders in opacity with full geome- 
try will require computational power beyond what is cur- 
rently available, we first present an analytic proof that a 
cone-like pattern of medium-induced emission persists to 
all orders in the correlations between multiple scattering 
centers (elementary emitters) and we focus on the case of 
immediate interest: light quark and gluon jets and final- 
state bremsstrahlung. Generalization to massive partons 
can easily be achieved, see e.g. [l^, but it is important 
to note that the effect of a heavy quark mass versus the 
jet energy depends on the coherent scattering regime 
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A. Radiative energy loss in the GLV formalism 

In our calculation we will use the GLV formalism of 
expanding the medium-induced radiation in the correla- 
tions between multiple scattering centers @. We first 
recall the definitions of the Hard, Gunion-Bertsch and 
Cascade propagators in terms of the gluon transverse mo- 
mentum k and the transverse momentum transfers from 
the medium q^: 



H 

B, 



' ^{ili2---i 
H — Ci , B(,j 



(k - qi^ - qi 



) 



(k 

^)Ul■■ 



(27) 



The relevant inverse gluon formation times can be writ- 
ten as: 



-'(iii2---im) 



= [k- 



(28) 



For final-state radiation, the intensity spectrum reads: 



^ dNs{FS) ^ CrOs 
dk+d^k 7r2 



E 



n 



dAzi 



n / d% 



1 daeliZj) 2 



(Telizj) d^qj 



-2C(i_..._„) • ^ B(„,+i_...^„)(„_...^„) cos ^W(fc^...^„)Azfc -cos ^ajf^k.-- ,n)^Zk 



,(29) 



, k=2 



, k=l 



where ^2 = ^^'^ '^{n+i,n) = B„ are understood. In 
the case of final-state interactions, zq w is the point of 
the initial hard scatter and = L is the extent of the 
medium. The path ordering of the interaction points, 
zl > 2j+i > Zj > zq, leads to the constraint — 
zl ■ One implementation of this condition would be Azi G 
[0, zl — ^-^j ] ^1^*^ it is implicit in Eq. (j29p . 

There is an obvious limit of the GLV radiative spec- 
trum when L 3> Ag 3> Tj, where Xg is the mean free 
path of the gluon in a hot QGP. Here, the contributions 
of the cos(- • • ) terms vanish after integration over the un- 
observed qi or Azi due to rapid oscillation. It is easy to 
see in this limit for n=l that. 

In the very high energy limit _B — s- 00, leading to large 
k phase space, a change of variables k — > k — qi shows 
that the first two terms in Eq. (|5D|) . cancel, leading to 
an incoherent Bertsch-Gunion gluon emission in a hot 
QGP medium with (n) — By direct inspection one 
can see that the n > 2 terms do not contribute. In fact, 
it is easy to verify that for any bremsstrahlung regime, 
initial-state, final-state and no hard scattering, this limit 



holds 0. More generally, in this limit it can be shown 
that the reaction operator R ^ 0. Naturally, for finite jet 
energies there will be corrections when k~^dN^ /dk^d'^'k is 
evaluated numerically with actual kinematic bounds 0]. 



B. Collinear radiation in GLV formalism 

While the example given above illustrates that limits 
can be imposed and taken in the GLV results, such limits 
are artificial in that the formation time of the gluon at 
the emission vertex spans Tf £ (0, 00). The Reaction 
Operator approach [61], i.e. the GLV formalism, is not 
an approach of averages: it compares differentially Tf 
to the separation between the scattering centers. For 
example, even when k — > the formation time can be 
small or large, depending on the momentum transfers for 
the medium. Let us investigate this case in more detail: 
we note that A;"*" « 2lu and k « rwn, where r is the 
angle relative to the jet axis. Here, n is a unit vector 
transverse to the jet axis which defines the azimuthal 
angle (f> of gluon emission. Using the results of Eq. (|29p . 
the 2D {(f>, r) angular distribution of gluons at n-th order 
in the correlated scattering expansion reads: 
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LodNS 

lim — ; — — CX UJ 

r^o dujd(j)dr 



1 daei{zj) 



i{zj) d'^q. 



qi 



(qi 



E 



qm+i 



q,i 



(q™ H h q„)2 



X cos 



E 



(qfe + 



2tj 



-Azk ] — cos j 

/ \fe=i 



(qm+i H h q,i)^ 

(qfc + • • • + q„)' ^ 



2u; 



-Azk 



(31) 



Here, we have already set r — where possible. We 
can use this general notation as long as we clarify cer- 
tain special cases: for to = n we have cos[(q„+i + 
q„)^Az„_|_i/2a;] = 1. For the transverse propaga- 
tors and TO we have a;r(q„+i -I- qn)/(qn+i -I- qn)^ = 
n. It is know that the leading n = 1 contribu- 
tion to final-state medium-induced radiation leads to 
limr^Q ojdN^ /dujd(j)dr — [22]. Our goal is to show that 
this result is general and holds to any order in the expan- 
sion. Its implications are that there is very little overlap 
between the techniques used to compute the "vacuum" 
and medium-induced contributions to the jet shape. A 
general proof requires demonstration of the absence of 
unprotected divergences for any set of momentum trans- 
fers {qi}, finiteness of the momentum transfer integrals 
as qi ^ oo and a mechanism that kills the small-angle 
contribution. 

1. We first look at the large qi limit. The transverse 
propagator contribution itself in Eq. ([3T|) behaves 
as ^ l/qf. Furthermore, irrespective of the small 
qi behavior of the momentum transfer distribution 
from the medium, for large momentum transfers 
the collisional cross section is suppressed by the 
Rutherford ~ 1/qf behavior, ensuring the finite- 
ness of the integrals. 

2. Next, we examine the potential singularity as |qi + 
■ ■ ■ + qn| ^0. The difference in the LPM interfer- 
ence terms in this limit goes as ©((qi H- • • • -I- q„)^) 
and for the most problematic transverse propaga- 
tor term (m=l) even as 0((qi + • • • + qn)"*). In 
summary, not only is there no divergence in this 
case, but the integrand in Eq. (I3T|) vanishes. 

3. We now collect the interference phases associated 
with problematic propagators as |qA: + - • •+qn| — ^ 0, 
1 < k < n. Expanding for a small net transverse 
momentum sums we find that the singularity is can- 
celed: 



fc-i 
si^E 

k-l 

-sin^ 



(qj 



qn)' 



2uj 



(qj+--- + q„)- 
2lj 



Az 



(qfe 



qn)' 



2lu 



(32) 



Azfc 



Actually, the lack of singularities persists also away 
from the small r limit. 



4. With the integrand well behaved and all integrals 
finite we see that the phase space factor r in the 
numerator is sufficient to ensure vanishing medium- 
induced bremsstrahlung contribution at the center 
of the jet. It is assisted by partial cancellation from 
angular integrals of the type / qi • q^ /(qi, qj) dcpij. 
It is only for the special case of h ■ (qi + • • • + qn) 
where the antisymmetric integrand under qi 
— qi for all i fully ensures the vanishing zero-angle 
radiative contribution. 



This completes our proof that at any order in opacity 
ujdN. 



lim 



mod 



+0 divdcftdr 







(33) 



Numerical simulations, using Monte-Carlo techniques, 
confirm independently that dl^ /duid^k vanishes as k ^ 



C. Numerical methods and QGP properties 

Results relevant to the LHC phenomenology are cal- 
culated using full numerical evaluation of the medium- 
induced contribution to the observed jet shapes and the 
modification of the in-medium jet cross sections. Jet pro- 
duction, being rare in that (j(Et > Et min)?Ayi(&) ^ 1, 
follows binary collision scaling ^ (P N\,inJ d^yL±. In con- 
trast, the medium is distributed according to the num- 
ber of participants density ^ d^TVpart./c'^xi^. Soft par- 
ticles that carry practically all of the energy deposited 
in the fire ball of a heavy ion collision cannot deviate a 
lot from such scaling. We take into account longitudi- 
nal Bjorken expansion since transverse expansion leads 
to noticeable corrections only in the extreme (3t ^ I 
limit 24]. In our approach all relevant finite time and fi- 
nite kinematics integrals, such as the ones over the sepa- 
ration between the scattering centers Azi = Zi — Zi^i, the 
bremsstrahlung gluon phase space Aqcd < < -Ejot., 
^QCD < k± < 2ll! 43], and the transverse momentum 
transfers < gi < ^s/4 — ^TO£)i?jet/2, are done nu- 
merically [H]. In our simulation we generated in-plane 
jets, (/)jet - (/'reaction plane = 0. This is of little impor- 
tance in central Pb-|-Pb collisions (b=3 fm), where the 
medium effects on jet propagation are most pronounced, 
but in semi-central (b=8 fm) and peripheral (b=13 fm) 
reactions this will lead to smaller than average energy 
loss. 
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The evolving intrinsic momentum and length scales in 
the QGP expected to be created at the LHC are deter- 
mined as follows: we first estimate the QGP formation 
time To = 1/(pt) = 0.23 fm, where (pt) ~ 850 MeV 
is obtained from extrapolations to LHC energies made 
by using Monte Carlo event generator results, fit to the 
CDF collaboration data from =1.8 GeV p + p col- 
lisions Here we account for the observed ~ 25% 
increase in the mean transverse momenta in going from 
N + N to A + A collisions at RHIC. Gluons domi- 
nate the soft parton multiplicities at the LHC and their 
time- and position-dependent density can be related to 
charged hadron rapidity density in the Bjorken expansion 
model (H: 



_ 1 (f{dN9/dy) 



1 3 
r 2 



drj 
dy 



d^idN'^^/dr]) 



(34) 



Here, dN'^^/d-q = KiVpart./2 with k w 9 for = 5.5 TeV. 

Table H] summarizes characteristics of Pb+Pb collisions 
at the LHC and initial QGP properties. An inelastic 
cross section (Ti„ = 65 mb has been used in an optical 
Glauber model where necessary. Assuming local thermal 
equilibrium one finds: 



T{t, xx) - V^'p(^,xx)/16C(3) , r > To 



(35) 



The Debye screening scale is given by rriD = gT, re- 
calling that we work in the approximation of a gluon- 
dominated plasma and Nf — 0. The relevant gluon 
mean free path is easily evaluated: Ag = l/a^^p with 
a99 ^ (9/2)7ra^/m|,. Note that in Table |T] the quoted 
initial mean temperatures, Debye screening scales and 
gluon mean free paths are obtained as averages with the 
binary collisions weight Taa (x ; 6) . In our evaluation 
we use gs — 2.5, (a^ — 0.5) to describe the scatter- 
ing of the jet with the medium but the QGP-induced 
bremsstrahlung is calculated with a running as{kT) for 
emission vertex, similar to the MLLA approach for the 
vacuum jet shapes. 

Evaluation of the medium-induced energy loss and its 
contribution to jet shapes is numerically expensive. Ex- 
act results have been obtained only for i^jot = 20, 100 
and 500 GeV. We interpolate for other values of interest. 



Class 


Central 


Mid-central 


Peripheral 


h[fm] 


3 


8 


13 


N part 


361 


165 


18 


dNO /dy 


2800 


1278 


137 


{T{ro)) [MeV] 


751 


693 


426 


(mo (to)) [GeV] 


1.89 


1.73 


1.07 


(A,(r„)) [fm] 


0.25 


0.27 


0.46 



TABLE L Summary of the relevant energy loss parameters 
ind initial QGP properties for central, semi-central and pe- 
ripheral collisions at -^s = 5.5 TeV collisions at the LHC. 



Jet axis 


R{2) 









/ / 


(i)""(2) 

















FIG. 4: (Color online) Schematic illustration of the cone 
radius R and the particle/tower pT / Et selection. The 
measured energy is the one that comes from particles with 
Pt > uimin and within R. 



D. Energy loss distribution 

A consistent energy loss theory provides complete in- 
formation for the differential distribution of the lost 
Ai^rad, i-e. the bremsstrahlung spectrum in Eq. (|29p . 
The main point that we make here is that this distribu- 
tion is completely determined by the properties of the 
QGP and the mechanisms of energetic quark and gluon 
stopping in hot and dense matter. Therefore, selecting 
different jet radii R and px min of the particles will signifi- 
cantly alter both the jet shape and the amount of energy 
lost by the hard parton which can be recovered in the 
experimental measurement. In contrast, we have seen in 
Fig. [2] that the jet shapes scale approximately as a func- 
tion of r/i?, i.e. they are independent of the selection 
of cone opening angle R. The jet cross section weakly 
depends on R, unless i? ^ 0. Finally, Zmin = 0.1 — 0.2 is 
necessary to noticeably alter the jet shape, implying that 
~ 20% of the parent parton energy has to be missed via 
Pt mill cuts to observe significant effects on V'-uac(?')- 

Experimentally, a clear strategy will be to use the 
leverage arms provided by R {= ij™"^ in the evaluation 
of the AE'rad) andpT min (= w""" in the evaluation of the 
AiSrad) to determine the distribution of the lost energy. 
This is illustrated schematically in Fig. [4] Theoretically, 
the first quantity to be calculated is: 



E 



iR" 



1 

E 



dio 



dr- 



dlB 
dujdr 



(36) 

We present in Fig.[5]this fractional energy loss for a quark 
jet and a gluon jet of energy Ej^t = 20 GeV inside a jet 
cone of radius and with acceptance cut w™'". In- 

creasing R"^^^ or decreasing w™'" we will recover more of 
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FIG. 5: (Color online) 3D plot for the ratio of the energy that 
a partons loses inside a jet cone of opening angle R"^'^^ with 
ui > to the total parton energy. We have chosen a jet of 
energy Ejet 

and varied the jet radius i^™™ and the acceptance cut oj™'". 
The upper surface is for a gluon jet and the lower surface is 
for a quark jet. 



0.2 0.4 0.6 0.8 1 

pjmax 



FIG. 6: (Color online) 2D projections of Eq. dMl The left 
panel shows the fractional energy loss dependence on the jet 
radius J?™^^ (oj™'" — 0) and the right panel shows this depen- 
dence versus the acceptance cut w™'" (i?""" = R°°). A gluon 
jet of Ejet = 20 GeV in fo = 3, 8, 13 fm Pb + Ph collisions at 
LHC was used as an example. 



20 GeV in 6 = 3 fm Pb+Pb collisions at LHC due to fragmentation): 



and/or u^' 



xplcading parton 

£j) — Hj. 



(38) 



r>h^ 
^AA ■ 



the parent parton energy, lost via gluon bremsstrahlung. 
We note that in Fig. [5] the mean energy loss was cal- 
culated as an average over the probability distribution 
P{e]E) [131, reflective of multi-gluon fluctuations: 



(37) 



For large fractional energy losses, such as in the illustra- 
tive example of a 20 GeV jet in central Pb + Pb collisions 
at the LHC, AEg/ AEq is much smaller than the asymp- 
totic rati o Ca ICp = 9/4 due to the kinematic constraint 
AE <E^. 

The separate dependence of A£'™(i?'"^'', w""")/£' on 
the cone radius and the momentum acceptance cut is 
more clearly illustrated in Fig. [6l We show central, 
mid-central and peripheral collisions, impact parame- 
ters & = 3, 8, 13 fm, respectively, in Pb + Pb reac- 
tions at LHC at nominal ^/s. We notice that, not sur- 
prisingly, the ratio AE™ {R"^^^ , uj"™^) / E goes down at 
larger impact parameters because the energy loss of the 
jet decreases in peripheral collisions. More importantly, 
at each impact parameter there is a variation of the 
amount of the bremsstrahlung energy, recovered in the 
cone. This is precisely the variation that will map on the 
^''aa(^'^^^^ w'"'") observable. For example, in the limit of 
a very small opening angle and/or large momentum cut 
to eliminate the QGP-induced radiation the suppression 
should approximate that of leading hadrons (up to dif- 
ferences arising from the possibly softer particle spectra 



One can see that the typical choices, R ~ 0.4 and w™™ — 
2 GeV are a good starting point to explore the variable 
quenching of jets. 



IV. TOMOGRAPHY OF JETS IN HEAVY ION 
COLLISIONS 

The purpose of this Section is to relate the theory of jet 
propagation in the QGP to experimentally measurable 
quantities. 



A. Experimental observables 

An essential ingredient that controls the rela- 
tive contribution of iPvu-cXt / R) and V'mcd.(?"/-R) — 
{l/AErad)dI^/dr to the observed differential jet shape 
in heavy ion reactions and also determines the attenua- 
tion of the jet cross sections is: 



/ ^ / 



R 



R 



AE, 



ad{(0,i?); 



mm _ Q 

,E)} 



A£;,.ad{(0,i?°°);(0,i;)} 



(39) 



the fraction of the lost energy that falls within the jet 
cone, r < R, and carried by gluons of a; > w™'" relative to 
the total parton energy loss without the above kinematic 
constraints. If this fraction is known together with the 
probability distribution P(e) for the parton energy loss 
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the medium-modified jet cross section per binary N - 
scattering can be calculated as follows: 



■N 



AA 



de 



e=0 



9:9 



1 



(1 - (1 - • e) 



d^E^dy 



(40) 



where E"^ 



Er/il - (1 - Ug) ■ e). The (1 - Ug) ■ e fac- 
tor accounts for the total "missed" energy in a jet cone 
measurement, which necessitates E^ > Et, and the Ja- 
cobian J — \d^ Elp / d"^ Et\ is properly accounted for. In 
this paper possible fluctuations of /g,g independent of e 
are not considered. Simple analytic limits illustrate the 
physics represented by Eq. ([40]) : if there is no energy loss, 
P(e) = i5(e), J = 1 and the cross section is unaltered. In 
the opposite limit, P(e) = 5{e — 1), the quenching of 
jets, if any, is completely determined by the fraction of 
the lost energy fq^g that is recovered in the experimen- 
tal acceptance. When fq^g = 1 once again there will be 
no attenuation of jets and when fq g — + our result ap- 
proximates the inclusive particle Raa{pt) [13]) see also 
Eq. (IMD. 

Next, we obtain the full jet shape, including the con- 
tributions from the vacuum and the medium-induced 
bremsstrahlung: 
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FIG. 



Aot.{r/R) = 



1 



Norm 



Je=0 



(Color online) Nuclear modification factor 
^^^-11 ,aj™'") as a function of the jet transverse energy, 
Et, at impact parameters 6 = 3 fm (circle), & = 8 fm 
(square) and 6 = 13 fm (diamond) in Pb-|-Pb collisions with 
^ = 5.5 TeV. 



same sum rule will hold in the presence of a medium 
since the total momentum of the partons must be con- 
served irrespective of whether they are propagating in 
vacuum or in the QGP with/without medium-induced 
bremsstrahlung (2^ . [27| . For jets, taking a monochro- 
matic pulse in the vacuum. 



d^E'dy 



(1 - (1 - h,9) 
[(1 - e) 4,^, (r/R) 



1 a 



NN 



— (5^(Et — Eq) , 



ad^Er 

we can easily verify that in the presence of a QGP 



mcd. 

We recall that, by definition, the area under any differ- 
ential jet shape, il)tot. {r/R), ipvac. {r/R) and V'mod. {r/R), 
is normalized to unity. Integrating over r in Eq. (j4ip . it 
is easy to see that the correct "Norm" is the quenched 
cross section, Eq. PO]) . The interested reader can inde- 
pendently carry out the analysis of the simple limiting 
cases and gain insight into the dominant contribution to 
the full jet shape. Proper treatment of isospin is implicit 
in Eqs. (gD]) and ((iT|l . 



B. Energy sum rule 

Sum rules provide useful integral representation of con- 
servations laws, originating from symmetries in QCD. 
One such example is momentum conservation in inde- 
pendent fragmentation: 



(41) 



d-'Ej 



1 cr^^(i?^ 00,0;" 



0) 



d'^Ej 







(43) 



(44) 



^/ zDh/q,g{z,Q'^)dz ^ 1 , 



(42) 



in case of perfect experimental acceptance. More gener- 
ally, only a fraction, 1 — (1 — /)(e), of Eq is recovered. 



C. Numerical results 

Combining our full theoretical model for the jet shape 
and the jet cross section in heavy ion collisions with re- 
alistic numerical simulations of parton propagation in 
the QGP, see Section IIII CI we first evaluate the nu- 
clear modification factor i?J^^(i?'"^'^, w™'") in Pb-f Pb col- 
lisions with center of mass energy = 5.5 TeV at 
the LHC. Fig. [7] illustrates the attenuation of the mea- 
sured jet rate as a function of the jet energy Et for 
different centrality classes. We use impact parameters 
b — 3, 8, 13 fm in conjunction with a jet cone radius 
^max ^ Q 4 g^^^ acceptance cut (w""" = GeV). The 

evolution of i?^^^(i?™'*'^, w™'") for jets is similar to the 



-kN. 



2/3 



part. 



where 



PT,h/P' 



'T,q{g) 



is the momentum fraction of 



parent partons carried by fragmentation hadrons. The 



one for leading particles in that In Raa 
However, k will depend on the selection of and w™'" 
in addition to the steepness of the underlying jet spectra 
and the properties of the QGP. We recall that the energy 
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FIG. 8: (Color online) iST-dependent nuclear modification 
factor i?j4*^(i?'"="',tj'"'") for different jet cone radii i?""''" (top 
panel) and at different acceptance cuts a;™'" (bottom panel) 
in 6 = 3 fm Pb+Pb collisions at ^ = 5.5 TeV. 



of the jet that will be redistributed out of the cone is 
e(l — /), see Eq. (|^. and the variation of its quenching 
with centrality is related to the fractional parton energy 
loss e = AE/E and its fluctuations, given by P{e)- 

Fig.lHldemonstrates the sensitivity of i?j^^(i?™'''', w™'") 
to the properties of the medium-induced gluon radiation 
through the independent variation of a;™'" and R^^^ ^ ad- 
vocated in this paper. For a fixed impact parameter, 
6 = 3 fm, the top panel shows a study of the quenching 
strength versus the jet cone radius when w™'" = GeV. 
In the approximations that we employ R™^^ = 2 is 
the upper bound of the medium-induced bremsstrahlung 
opening angle relative to the jet and, consequently, con- 
stitutes perfect experimental acceptance. In this case 
there is no deviation from binary collisions scaling. The 
smooth evolution of B}^j^{R"^'^^ , oj™™) with decreasing 
jjmax ^ signature of the large- angle gluon radiation pat- 
tern in the QGP Note that if dl^ /dojdr were pre- 
dominantly coUinear, there would be no deviation from 



For a finite R^"^"" = 0.7 even if w"'" = GeV 
j^^{R™^^, w™™) does not reach unity, see our discussion 
The largest variation in the qu' 

,min 



unity. For R™^^ < 0.2 the magnitude of jet quenching 
approaches the suppression for leading hadrons. A good 
starting point is a cone radius selection R'^^^ — 0.4 — 0.7 
if the experimental statistics allows for positive identifi- 
cation of 30% to a factor of 2 variation in the measured 
cross section. In the bottom panel of Fig. [5] we present 
the sensitivity of jet attenuation to the minimum parti- 
cle momentum/calorimeter tower energy deposition cut 

LU 

pj^t f pmax , ,min 

above. The largest variation in the quenching strength is 
observed between w™'"^ = 2 GeV and w™'" = 5-10 GeV, 
and reflects the typical energy of the stimulated gluon 
emissions. We emphasize that, in the GLV approach 
partons lose energy through few GeV bremsstrahlung 
gluons [27]. For w"'" > 10 GeV, w""") ap- 

proaches again the characteristic leading particle sup- 
pression. In summary, for the same centrality, Et and 
^/s the continuous variation of quenching values may help 
differentiate between competing models of parton energy 
loss , thereby eliminating the order of magnitude un- 
certainty in the extraction of the QGP density. 

Detailed investigation of i?^^^(i?™''^, cj™'") can also in- 
dicate whether "elastic" 2^2 processes, such as colli- 
sional energy loss [2§|, or "inelastic" 2 — > 2-|-n processes, 
such as bremsstrahlung 0] and hadron dissociation in the 
QGP [13, dominate the inclusive particle and particle 
correlation quenching. If the energy loss per interaction 
in the first scenario AEcoi\./E < 5%, the recoil parton 
form the medium will accelerate almost transversely rela- 
tive to the jet axis and will not be part of the jet for any 
reasonable selection of R"^'^^. Therefore, for coUisional 
energy loss, in contrast to the well-defined evolution of 
the jet suppression with cone radius and the acceptance 
cut seen in Fig. [51 the cross section attenuation will be 
large and constant and will approximate the quenching 
of leading hadrons. Note that R'Xli < 1 has also been ob- 
served in Monte-Carlo silumatios of jet quenching [sij. 

We now turn to the numerical results for the jet shape 
in Pb-t-Pb coUisions at ^ = 5.5 TcV at the LHC. In 
Fig. [5] we first explore the difference between the vac- 
uum and the medium-induced only {Et given for the 
parent parton) 'ip{r / R) as a function of the impact pa- 
rameter, jet energy, and the cone radius. We note that 
in central heavy ion reactions for lower Et and, in par- 
ticular, for R > 0.7 the two differential shapes can be 
quite similar. The differences become more pronounced 
for smaller jet radii where the experimental acceptance 
will subtend the part of phase space with the most effec- 
tive cancellation of the coUinear medium-induced radia- 
tion [2^ . It is interesting to observe that in going to more 
peripheral collisions tpmcd.(r / becomes slightly wider. 
The underlying reason is that the LPM destructive in- 
terference between the radiation induced by the large 
scattering and the radiation induced by the subsequent 
interactions in the QGP determines the angular distri- 
bution in the bremsstrahlung spectrum. Thus, a small 
medium size facilitates the resulting cancellation for glu- 
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FIG. 9: (Color online) The vacuum and medium-induced only 
jet shape for Et = 20 GeV (top panel) and Et = 100 GeV 
(bottom panel) at impact parameters b = 3, 8, 13 fm in 
Pb+Pb collisions at the LHC. Two jet cone radii, R — 0.7 
and 0.4 are shown. 



ons of large formation time. A general observation is that 
the medium tends to redistribute the flow of energy more 
evenly inside the jet cone, especially for large r/R 1 

The pattern of energy flow for in-medium jets is shown 
in Fig. [TO] together with ipvac.{r/R) and tpmcd. {r/R) for 
comparison. We used Et — 20 GeV to 200 GeV and 
R — 0.4 to 0.7 to cover a wide range of measurements 
that will become accessible during the first year of heavy 
ion running at the LHC. One observes that there is no 
significant distinction between the jet shape in the vac- 
uum and the total in- medium ip(r/R). The underlining 
reason for this surprising result is that although medium- 
induced gluon radiation produces a broader i'mcd.i''' / R) , 
this effect is offset by the fact that the jets lose a finite 
amount of their energy, see Figs. [5] and [51 Furthermore, 
when part of the lost energy is missed due to finite exper- 
imental acceptance, the required higher initial virtuality 
jets are inherently narrower, see Figs. [Hand [3) 

In Table HI] we show the mean relative jet radii (r/R) in 
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FIG. 10: (Color online) Comparisons of the jet shape in vac- 
uum, the medium-induced jet shape, and the total jet shape 
for cone radii R = 0.7 and R = 0.4 and four different energies 
Et = 20, 50, 100, 200 GeV, respectively, in central Pb-|-Pb 
coUisions at the LHC. 



the vacuum and in the QGP medium created at the LHC 
for two different cone selections R = 0.4 and R = 0.7 and 
four transverse energies Et = 20, 50, 100, 200 GeV. We 
see that in the realistic numerical simulation there is very 
little < 10% increase in the magnitude of this observable. 
The difference is slightly larger for a smaller cone, since 
it emphasizes the large-angle character of the medium- 
induced radiation [6','25|. Therefore, a rough 1-parameter 
characterization of energy flow in jets will not resolve the 
effect of the QGP medium. It can, however, exclude sim- 
plistic scenarios of full jet stopping in the QGP that lead 
to {r/R) growth by as much as 60%. It is also impor- 
tant to stress that the QGP is rather "gray" than "black" 
and only a fraction of the energy of the parent paron is 
lost via stimulated gluon emission. The effect of even a 
moderate A£;"'(i?™'^'', a;"™)/£; can be amplified by the 
steeply faUing cross sections for the i?^^^(i?T; R"^^^, w™'") 
observable, see Fig. [SI but this is not the case for {r/R). 

Lastly, we point out where the anticipated jet broad- 
ening effects will be observed in the differential shape by 
studying the ratio Aot. {r / R) /i/j^e.dr/ R) in Fig. [Ill We 
have used the same transverse energies and cone radii as 
in Fig.[Tni We recall that the small r/R < 0.25 region of 
the intra-jet energy ffow in p+p collisions in our calcula- 
tion has uncertainties associated with the normalization 
of the jet shape. In the moderate and large r/R > 0.25 
region our theoretical model gives excellent descriptions 
of the Fermilab Run II (CDF II) data, as shown in Fig.[Tl 
The QGP effects are manifest in the "tails" of the energy 
flow distribution and for a cone radius R = OA the ra- 
tio could reaches ~ 1.75 when r/R —f 1. However, for 
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R = OA 


Vacuum 


Complete E-loss 


Realistic case 


(r/R), Et = 20GeV 


0.41 


0.55 


0.45 


(r/R), Et = 50GeV 


0.35 


0.48 


0.38 


(r/R), Et = lOOGeV 


0.28 


0.44 


0.32 


(r/R), Et = 200GeV 


0.25 


0.40 


0.28 


R = 0.7 


Vacuum 


Complete E-loss 


Realistic case 


{r/R}, Et = 20GeV 


0.41 


0.44 


0.42 


(r/R), Et = 50GeV 


0.33 


0.39 


0.37 


(r/R), Et = lOOGeV 


0.27 


0.34 


0.29 


(r/R), Et = 200GeV 


0.24 


0.31 


0.26 



TABLE II; Summary of mean relative jet radii (r/R) in the vacuum, with complete energy loss, and in the QGP medium. 
Shown are results for cone radii R = 0.4 and R = 0.7 and transverse energies Et ~ 20, 50, 100, 200 GeV at ^/s — 5.5 TeV 
central Pb+Pb collisions at the LHC. 



experiments to observe this enhancement of the ratio of 



the total jet shape in medium to jet shape in vacuum at 
r/R > 0.5, high statistics measurements will be needed. 
This precision can hopefully be achieved with the large 
acceptance experiments at the LHC. 
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FIG. 11: (Color online) The ratios of total jet shape in heavy- 
ion collisions to the jet shape in the vacuum for jet energies 
Et = 20, 50, 100, 200 GeV. Two cone radu R = 0.7 (top 
panel) and R = 0.4 (bottom panel) at 6 = 3 fm in Pb+Pb 
collision with -v/s = 5.5 TeV were chosen. 



V. CONCLUSIONS 

The unprecedentedly high center of mass energies at 
the LHC will usher in a new era of precision many-body 
QCD. The theory and phenomenology of jets in nuclear 
collisions are expected to evolve as the new frontier in 
the perturbative studies of parton propagation in the 
QGP [2l[. In this paper we discussed three important 
aspects of such studies: a generalization of the analytic 
approach for calculating differential jet shapes [T^ that 
can accommodate experimental acceptance cuts needed 
to isolate jets in the high multiplicity environment of 
heavy ion collisions; the theory of the intra-jet energy 
flow redistribution through large-angle medium-induced 
gluon bremsstrahlung; and a comprehensive new set of 
experimental observables that can help identify and char- 
acterize the mechanisms of parton interaction in nuclear 
matter. 

In elementary nucleon-nucleon collisions we compared 
our theoretical model to the CDF II Tevatron data on 
jet shapes [l^l and investigated the basehne i^vacXr / R) 
at the higher = 5.5 TeV at the LHC. We found that 
in the absence of a hot and dense QGP matter these 
shapes are self-similar and approximately independent of 
the cone radius R. Elimination of low momentum parti- 
cles of up to '^few GeV is not likely to significantly alter 
the pattern of intra-jet energy flow for Et > 50 GeV jets. 
In nucleus-nucleus reactions we demonstrated that the 
characteristic large-angle QGP-stimulated gluon emis- 
sion [l^l persists to all orders in the correlation between 
the elementary bremsstrahlung sources. We showed that 
this intensity spectrum can be fully characterized by the 
amount of the lost energy that falls inside the jet cone 
(r < i?™'''^, u! > w™'") and derived the medium modifica- 
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tion of the jet shapes and jet cross sections in the QGP, 
subject to an intuitive energy sum rule. 

To demonstrate the connection between the QGP 
properties, the mechanisms of parton interaction and en- 
ergy loss in hot and dense matter, and a new class of jet- 
related experimental observables, we carried out realistic 
simulations of quark and gluon production and propaga- 
tion in the medium created in relativistic heavy ion colli- 
sions at the LHC. We introduced a natural generalization 
of the leading particle suppression to jets and showed 
that it is a more differential and powerful tool that can 
be used to assess in approximately model-independent 
way the characteristic properties of the induced gluon 
intensify spectrum. Consequently, in the future progress 
can be made toward identifying the set of approxima- 
tions [1, H, 0, 0] that most adequately reflect the dy- 
namics of hard probes in the QGP. We also discussed 
how the evolution of R^^KR'^^'', w"'") with the jet cone 
radius R™^^ and the acceptance cut w™™, or the lack 
thereof, can help differentiate between radiative and col- 
lisional energy loss paradigms of light and heavy quark 
attenuation. The theoretical approach, developed in this 
manuscript, allows to investigate the correlation between 
the quenching of jets and the in-medium modification of 
their shape. Surprisingly, up to five-fold attenuation of 
the cross section corresponds to a rather modest < 10% 
growth of the mean relative cone radius (r/R). The an- 
ticipated broadening of jets is most readily manifest in 
the periphery of the cone, r/R 1, and for smaller radii, 
e.g. i?'"^'^ = 0.4. 

Further refinements in jet phenomenology, especially 
the consideration of jet cross sections, should include cold 
nuclear matter effects, such as nuclear shadowing, the 
Cronin effect, and initial state energy loss [s^]. We finally 
note that the study of inclusive jet shapes and cross sec- 
tions in heavy ion collisions can easily be generalized to 
hadron, photon or di-lepton tagged jets O [H, H^l with 
the benefit of additional constraints on the hard process 
virtuality and the parton energy. 
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APPENDIX A: JET CROSS SECTIONS 

In this paper we focus exclusively on large momentum 
transfer processes, » ^qcD' that can be systemati- 
cally calculated in the framework of a reliable theory, the 
perturbative QCD factorization approach. Factorization 
not only separates the short- and long-distance QCD dy- 



namics but implies universality of the parton distribution 
functions (PDFs) and fragmentation functions (EEs) and 
infrared safety of the hard scattering cross sections. Eor 
hadronic collisions, one of the most inclusive processes 
is jet production. To lowest order (LO) the invariant 
differential cross section reads [35| : 

Vd max 

d<Th^ht _ j^y^ f <Pa/hS^a,^^f)(t>b/hS^b,^^f) 

dycCppT^ ^ J XaXb 

abcdy^ 

x^%^|M.,^.,r (Al) 

Here, s = (Pa + Pb)'^ is the squared center of mass energy 
of the hadronic collision and Xa — pt/Pai ^b = Pb l^b 
are the lightcone momentum fractions of the incoming 
partons. In this formulation, for masslcss initial-state 
quarks and gluons, 

y.max(min)=+(-)lnf^-^e+(-)^=) , (A2) 
\mTd rriTd J 

where m^. = mj +p^.. In Eq. (jAl[) 4>i/}^.{xi, i) is the 
distribution function of parton "i" in the hadron hi and 
/ir and /i f i are the renormalization and factorization and 
scales, respectively. In this work calculations are done 
strictly in the coUinear factorization approach and we use 
the CTEQ6.1 LO PDFs \Mab-^cd\'^ are the squared 
matrix elements for ab — s- cd partonic sub-processes. 

Numerical results for inclusive jet cross sections in high 
energy hadronic collisions are shown in Fig. 1121 here 
Pt = Et- The top panel compares the LO calculation, 
Eq. (jAl[) . to CDF data on jet cross sections in p -I- p at 
\/s = 1.96 TeV. Excellent agreement between data and 
the theory using K = 1.5, independently extracted from 
the charged hadron + differential cross section at 
the Tevatron. It indicates ^ 50% next-to-leading order 
correction. Alternatively, we have studied the sensitiv- 
ity of the cross section to the choice of the factorization 
and renormalization scales by varying fir = fJ-f = Pt/'^, 
Pt and 2pT. Not surprisingly, the uncertainty is also on 
the order of ~ 50%, similar to the phcnomenological K- 
factor. Jet shapes depend on the parton species, quarks 
versus gluons, and the insert shows the fraction of quark 
jets versus px at the Tevatron. The bottom panel gives 
predictions for the corresponding jet cross sections at the 
LHC per nucleon-nucleon collision for y/s = 14 TeV and 
5.5 TeV, without quenching. Insert shows the increased 
fraction of gluon jets relative to the Tevatron. 

We can now evaluate the feasibility of differential jet 
shape measurements at the LHC. During the first three 
years of running, even at a fraction of the designed 
C = 10'^'* cm~^s~^, LHC is expected to deliver an in- 
tegrated luminosity of 10 fb~^ per year. In heavy ion 
collisions, nominal £ = 10^^ cm^^s"^ is not expected to 
be achieved either. An integrated luminosity of 1 nb~^ 
per year is a realistic projection. As shown in this pa- 
per, the anticipated quenching factor for energetic jets 
depends on the selection of R™^^ and w™™. In the limit 
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FIG. 13: (Color online) Comparison of the theoretically com- 
puted jet shapes with Raep = 1-3, 2 and two different non- 
perturbative correction scales Qo = 2.0 GeV and 3.0 CeV to 
experimental data [20|l. A jet cone radius R = 0.7 was chosen 
in ^ = 1960 CeV p + p collisions by CDF II. 



FIG. 12: Top panel: Inclusive jet cross section in p+p col- 
lisions at the Tevatron ^/s = 1.96 GeV calculated to LO in 
PQCD and compared to the CDF run II data [33|- Insert 
shows the fraction of gluon jets as a function of pr- Bottom 
panel: predicted baseline jet cross sections in p-l-p collisions 
at the LHC at = 5.5 TeV and 14 TeV. 



of narrow jets R^^^ w i?^^ = 0.25 - 0.5 [13, HI]. Tak- 
ing this into account, but neglecting for the moment the 
complications associated with jet reconstruction in the 
high particle multiplicity environment of heavy ion col- 
lisions [l3], from Fig. [T^we find that excellent < 10% 
statistical precision can be achieved for inclusive mea- 
surements of jets of Pt as high as 160 GeV in Pb-I-Pb re- 
actions and 1.3 TeV in p-f p reactions. Jet shape measure- 
ments require higher statistics since the shape functions 
are precipitously falling as of r/i? —> 1. We expect that 
very good, < 30% at large r/R ~ 1, jet shape measure- 
ments will be possible to pr as high as 100 GeV and 900 
GeV in Pb-I-Pb and p-l-p collisions, respectively. These 
will give an indication as to whether medium-induced 
broadening is present in the tails of the intra-jet energy 
distribution. The estimates presented in this appendix 
are conservative. Ay = 1. 



APPENDIX B: CONTRIBUTIONS TO THE 
VACUUM JET SHAPE 



It is important that a jet finding algorithm be in- 
frared and collinear safe. In full Monte-Carlo simula- 
tions of high-energy hadronic events the algorithm can be 
"tested" and matched to the experimental measurement 
techniques. In analytic calculations an approximate way 
to mimic the effect of jet splitting/merger is to introduce 
an adjustable parameter, Rsep, for cone type algorithms. 
If two partons are within an angle IL^R of each other, 
they should be mer ged into one jet [3a |. This approach 
may not be optimal 1 1^ since it does not generalize intu- 
itively for NLO jet shape calculations. When comparing 
theoretical results to experimental data one finds that 
Rsep is a function of the event kinematics, i.e. it is jet 
momentum (energy) dependent. Nevertheless, at low- 
est order this is a useful phenomenological approach to 
obtain the best possible description of the baseline differ- 
ential jet shapes in nuclcon-nuclcon collisions, needed for 
the study of QGP-induced effects. Also, it known that 
for jet cross sections at NLO the results from other jet 
finding algorithms, such as the fully infrared and collinear 
safe kx algorithm with a jet-size parameter D, coincide 
within a few % with the results from the cone algorithm 
with appropriate matching/choice of R and Rsep ^16] , 
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FIG. 14: (Color online) Differential jet shapes in p + p col- 
lisions at y/s = 1960 GeV [S^l are compared to theory. 
Leading-order without initial-state radiation (LO w/o IS), 
leading-order contibution (LO), leading-order with power cor- 
rection (LO-I-PC), and leading-order with power correction 
and Sudokov resummation (LO-I-PC-I-RS) results are shown 
separately. 



In Fig. [13] we demonstrate how the numerical results 
for jet shapes can be optimized using experimental data 
at the Tevatron. Proper normalization of 7/;(r, R) in 
this Appendix is achieved via first bin subtraction and 
does not affect the moderate and large r/R part of 
the energy flow distribution. At low transverse energy, 
Et = 45 — 55 GeV, a numerical calculation with large 
Rsep = 2.0 gives a much better fit than one with a small 
Rsep = 1.3. At high Et = 250 - 277 GeV, the theoreti- 
cal result with small Rsep = 1.3 agrees with the data by 
CDF II fairly well. This finding is consistent with a pre- 
vious study showing that, to fit the data, at low jet Et 
one always needs a larger Rsep and that Rsep should drop 
with increasing transverse energy [38| . In the numerical 
calculations shown in Fig. [13] the contribution of power 
correction to the jet shape has been included. The scale 
Qo is used to separate the non-perturbative effect from 
the perturbative derivation (see Eqs. (|2ip and ([H])). In 
Fig. [13] two different values of the non-perturbative cor- 
rection scale, Qo — 2 GeV and Qo = 3 GeV, are used. 
It is clear that the curves with these two different scales 
are practically indistinguishable, which demonstrates the 
consistency of our treatment of non-perturbative effects. 

In Fig. [13] we illustrate the influence of the different 
perturbative and non-perturbative contributions to the 



jet shape for different Et- We observe that the effect 
of initial-state radiation (IS), absent in lepton colliders, 
is sizable. In fact, at very high jet energy, a leading- 
order (LO) calculation with initial- state radiation al- 
ready gives a good description of the experimental data. 
However, below « 75 GeV it is impossible to fit the 
data using only leading order, even with the maximum 
Rsep = 2. Other contributions should be considered in 
an improved theoretical description of jet shapes. These 
include the effect of the running coupling constant in the 
momentum transfer integrals (MLLA) and power correc- 
tions (PC) oc Qo/Et- Sudakov resummation (RS) en- 
sures the finiteness of tp{r, R) in the region r ^ 0. The 
explicit formulas for these two contributions are given in 
Section [il Al It can be seen from Fig. [TJ] that all pertur- 
bative and non-perturbative effects should be taken into 
account if reliable description of the experimental results 
is to be achieved. 



APPENDIX C: DOUBLE DIFFERENTIAL 
MEDIUM-INDUCED JET SHAPE 

In this paper we have investigated extensively the dif- 



ferential jet shapes ■ip{r, R) 



dr 



in vacuum and 



in heavy-ion collisions at the LHC. This quantity inte- 
grates over the energy distribution of the partonic jet 
fragments. Thus, the information about the angular dis- 
tribution of soft vs hard shower partons is lost. At the 
LHC it might be possible, via particle tracking or jet re- 
analysis, to recover this information on an event-by-event 



basis and construct 

C?^^ined.(?^, R) 

drdz 



drdz ' 



where 



1 dia{uj = zE. 



jet: ' 



AS"(i?,0) 



duidr 



(CI) 



In Eq. (|Cip the QGP-induced double differential shape is 
normalized such that it integrates to unity with R"^^^ — 
R, tj™'" = and z = u/Ejet- To illustrate the additional 
insight that can be gained through such studies we show 
numerical results for a simplified case where a quark jet of 
Ejet = 100 GeV propagating through a QGP of length 
L — 6 fm. While the medium is not expanding and 
characterized by Xg = 1.5 fm, itid = 0.7 GeV and as = 
0.3, the (AE/E) approximates the full numerical result 
for central Pb-fPb collisions at the LHC. 

Fig. [T5l shows a 3D plot of ty^^^^tt-^ for momentum 

^ ' ' ^ {27rrjdrdz 

fraction z = k+/E+ - LUgiuon/ET = 0.01, 0.03, 0.1, 0.3 
and i? = 1, w'"'" = GeV. We can observe that at 
z — 0.01 {oj — 1 GeV) the double differential medium- 
induced jet shape is dominated by gluon radiation at 
large opening angle r/R. At z = 0.03 {uj — 3 GeV) the 
peak lies in the intermediate r/R region with significantly 
suppressed gluon radiation at small open angle r/R. In- 
creasing Lu further narrows the medium-induced inten- 
sity profile. It is tempting to associate the characteristic 
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FIG. 15: (Color online) 3D plot of the double differential medium-induced jet shapes for a quark jet of jet energy Et = 100 GeV 
with R — 1, oj™'" = GeV in Pb+Pb collisions with ^/s = 5500 GeV at the LHG. Four different figures represent the double 
differential jet shape for medium-induced gluon momentum fraction z — 0.01, 0.03, 0.1, 0.3, respectively. 



shapes for small values of z in Fig. [15] with measurements 
of enhanced away-side large-angle particle-triggered cor- 
relations [39', 13, SH . However, the current RHIC data 
presents challenges in separating the jet from the back- 
ground or even distinguishing between events with 1 or 
2 minijets (recall that Ncoii ~ 1000 in central Au-|-Au). 
Analysis on an event-by-event basis can help reveal un- 



ambiguously the QGP medium response to jets. When 
future experimental measurements of jet in heavy ion 
collisions are perfected at the LHC full numerical sim- 
ulations of the double differential shape, including the 
vacuum and medium- induced components, will soon fol- 
low. 
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